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A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket- 
based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed 
within a statistical framework, commonly known as design of experiments. The parametric design space included 
four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A 
D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) 
solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear 
effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were 
executed with the NPARC v3.0 code. The response used in the statistical analysis was created from /sp efficiency 
data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both 
one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, 
iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant 
effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent 
that the free-expansion process, directly downstream of the rocket nozzle, can influence the / 5p efficiency. Numerical 
schlieren images and particle traces have been used to further understand the physical phenomena behind several 
of the statistically significant results. 
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Nomenclature 

area 

primary rocket nozzle area ratio (see Fig. 1) 
mixer-ejector inlet area ratio (see Fig. 1) 
mixer-ejector exit area ratio (see Fig. 1) 
constant coefficient of statistical model 
diameter 

vacuum specific impulse 

mixer-ejector nozzle length (see Fig. 1) 

mixer-ejector length/diameter ratio (see Fig. 1) 

reference length 

propellant mass flow 

secondary mass-flow injection ratio 

rocket chamber total pressure 

static pressure 

Reynolds number at throat conditions 
rocket chamber total temperature 
axial thrust 
velocity 

independent variable of statistical model (see 
Table 1) 

transformed independent variable [see Eq. (1)] 
grid refinement ratio 

spatial order of accuracy of computational fluid 
dynamics solution 
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t)f = vacuum specific impulse efficiency 

p = density 


Subscripts 


cav 

P 

r 

s 

3 

6 


annular cavity 
primary 
rocket nozzle 
secondary 
mixer-ejector inlet 
mixer-ejector exit 


Superscripts 

CFD = calculated from the computational fluid dynamics 
flow solution 

MAX = calculated from isentropic flow theory 
RSM = calculated from response surface model 

* = rocket throat 


I. Introduction 

R EUSABLE launch vehicles have many potential economic 
benefits. Currently, several single-stage-to-orbit propulsion 
concepts are under investigation within NASA’s Advanced Space 
Transportation program. One of the most promising technologies re- 
ceiving renewed attention is rocket-based combined-cycle (RBCC) 
propulsion. 

RBCC propulsion systems have been studied for many years. 
Much of the analysis to date can be traced back to the foundation es- 
tablished three decades ago. 1 Currently, several organizations, aca- 
demic and industrial, are investigating the possibilities of RBCC 
propulsion. 2-4 Olds and Lee 5 and Olds 6 have conducted qualitative 
multidisciplinary design optimization and economic analysis for a 
variety of RBCC vehicles. Often, such analyses will hinge upon 
assumed values of engine component efficiencies. However, a truly 
robust design optimization demands further investigation into the 
subsystem performance of these complex propulsion cycles. 
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The RBCC concept considered here is defined by four separate 
operational modes in a single-stage-to-orbit configuration. First, the 
engine functions as a rocket-driven ejector (mode 1). Then the rocket 
engine is switched off, and subsonic combustion is present in the 
ramjet (mode 2). As the vehicle continues to accelerate, supersonic 
combustion occurs in the scramjet (mode 3). The rocket is eventually 
reignited (mode 4) for the final ascent to orbit in an all-rocket mode. 
The overall performance of the combined-cycle engine is heavily 
dependent upon the efficiency of this all-rocket mode because a 
significant portion of the propellant mass is consumed during this 
final ascent mode. The present study focused upon the engine per- 
formance of this fourth and final mode. Engine performance was 
measured in terms of the thrust generated from a given amount of 
propellant, or vacuum specific impulse, / sp . This value was com- 
pared to the maximum achieved via isentropic expansion to obtain 
values of / sp efficiency. 

Mode 4 can be simulated with the familiar geometry of a ducted 
rocket. Geometries of this type have been studied extensively for 
applications concerning rocket driven diffuser flows, 7 base-pressure 
studies, 8 and supersonic sudden expansions. 9 However, the present 
study sought to quantify the expansion efficiency of a ducted rocket 
flow, exhausting to a vacuum, over a wide range of geometric and 
flow conditions. 

We examined the performance of RBCC systems operating in 
mode 4 with computational fluid dynamics (CFD). The primary 
reason for choosing numerical simulation over physical experimen- 
tation was the inherent cost advantage for a test program of this 
magnitude and complexity. NPARC v3.0 (Ref. 10) was chosen be- 
cause of the extensive validation for propulsion flows, including 
ducted rockets 11,12 and the variety of simulation options. 

We chose to execute the CFD simulations within the context of 
statistical Design of Experiments (DOE). 13 Problem complexity ne- 
cessitated this approach. The application of DOE methods within 
the aerospace research community is still relatively new, 14 despite 
being well established in other fields. Giuntaet al. 15 applied a novel 
DOE method to optimize the configuration of a supersonic trans- 
port aircraft using a hierarchy of aerodynamic models. Knill et al. 56 
applied response surface methodology (RSM) designs to optimize 
the configuration of a high-speed civil transport aircraft using Euler 
flow analysis. Tolle 17 applied RSM and Navier-Stokes flow analysis 
to optimize the definition of a hypersonic wind-tunnel nozzle con- 
tour. We applied RSM designs and Navier-Stokes calculations to 
analyze RBCC nozzle efficiency over a broad range of parameters. 

We varied six independent parameters simultaneously to study 
their effects upon system performance. A full parametric study of 
these six variables, each at three levels, would necessitate an enor- 
mous computational effort. DOE methods minimize the number of 
experiments, based upon the anticipated relationships between the 
response and the independent variables under investigation. This 
approach yielded a quantitative analysis of a generic RBCC nozzle 
performing in all-rocket mode. The significant factors influencing 
this performance have been estimated for future experimental anal- 
ysis. In addition, the present effort will help the propulsion system 
designer to focus upon the most important aspects when engaged in 
a trade study. 


II. Nozzle Configuration 

A generic RBCC nozzle operating in mode 4 can be represented 
by three simple components: the primary rocket engine, the annular 
secondary flow cavity, and the mixer-ejector duct (see Fig. 1). The 
relevant flow features can be represented by two propellant flow- 
paths: the primary flow through the rocket thrust chamber and the 
relatively small secondary flow injected into the annular cavity re- 
gion. The opportunity for utilizing a secondary mass-flow injection 
scheme would arise from the specification of a gas-generator engine 
cycle. For the present study a generic RBCC nozzle configuration 
was parameterized with the following two fiowfield and four geo- 
metric variables: rocket chamber total pressure P c , secondary mass- 
flow injection ratio mjm p , mixer-ejector inlet area ratio A 3 /A*, 
mixer-ejector length/diameter ratio L/D 3y mixer-ejector exit area 
ratio A 6 /A 3 , and primary rocket nozzle area ratio A r /A*. These 


Table 1 Independent variables: definition 
of the six -dimensional design space 


Independent 

variable 

Physical 

variable 

Low Mid High 

x\ 

P, (psi) 

300 750 

1200 

X2 

(m s jm p ) x 100% 

0 4 

8 

X3 

(A 3 /A*) 

40 120 

200 

X4 

(L/D 3 ) 

2 3.5 

5 

X5 

(A 6 /A 3 ) 

l 1.5 

2 

X6 

(A r /A*) 

4 12 

20 



Fig. 1 Axisymmetric RBCC nozzle geometry: A r denotes the rocket 
exit area, A 3 denotes the mixer-ejector inlet area, L denotes the mixer- 
ejector length, As denotes the mixer-ejector exit area, and A 41 (not shown) 
denotes the primary rocket throat area. 


assumptions led to the parametric design space depicted in Fig. 1 
and Table 1. We have intentionally chosen a broad design space for 
investigation so that the resulting statistical model of nozzle effi- 
ciency will encompass a wide range of performance variation. The 
transformation shown in Eq. ( 1 ) was used to scale each independent 
variable to a (— 1 , + 1 ) range. 


ixi - (*; ugh + *, iow ) 

(*, Weh -*, 1 ™) 


(i) 


Thus jtj e [xf m , jc)™ 1 , jc l h ' 8h ], whereas Jc, e [— 1 , 0, + 1 ]. This trans- 
formation simplified the statistical analyses presented later. 


III. Nozzle Efficiency 

The integrals for axial mass flow and vacuum thrust were cal- 
culated with a trapezoidal integration across the exit of the duct. 
The ratio of these integrals defined the specific impulse for a given 
simulation. 


m = /, 6 (p « 2 + p) 

f A Jpu)dA 


jC FD 

^ = 7 ^* 100% < 3 > 

sp 

Isentropic flow theory defined the maximum specific impulse that 
a given configuration can achieve. Thus the ratio of the calculated 
and isentropic impulse values defined the impulse efficiency. The 
isentropic analysis accounted for the flow through the primary' rocket 
engine alone; any additional injected flow is not reflected in the 
denominator of the efficiency. An alternate definition and analysis, 
based upon the gas generator cycle, was presented in Ref. 18. 

IV. Numerical Modeling 

A. NPARC v3.0 

The CFD simulations were performed with version 3 of the 
NPARC code configured to solve the axisymmetric Navier-Stokes 
equations. Two different turbulence models were examined so that 
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Table 2 Combustion chamber specifications for the 
three total pressure values 


Physical 

variable 

Condition 1 

Condition 2 

Condition 3 

P c psi 

300 

750 

1,200 

7o°R 

6,200 

6,400 

6,500 


321,000 

323,000 

324,000 


modeling sensitivity could be addressed. The two turbulence mod- 
els chosen for examination were the Spaiart-Allmaras (SA) one- 
equation model, 19 and the Chien ( k-E ) two-equation model for com- 
pressible flows. 20,21 

B. Boundary Conditions 

The thrust chamber inflow was specified with constant total pres- 
sure, total temperature, and flow angle. The physical walls of this 
nozzle were modeled with an adiabatic no-slip surface. The exit 
plane of the mixer-ejector was modeled with an extrapolation bound- 
ary condition to simulate the ambient vacuum environment. The 
nozzle axis was modeled with the symmetry condition. Finally, the 
secondary flow entering the annular portion of the duct was modeled 
with the fixed mass flow, total temperature, and flow angle. The sec- 
ondary mass flow entered the cavity with a uniform velocity profile 
applied across the annular upstream boundary. For cases executed 
without secondary flow, the fixed mass-flow boundary condition was 
replaced with the no-slip condition. 

C. Flowfield Parameters 

The reference length for these simulations was specified as the 
nozzle throat diameter (/ ref = 2.54 cm). The ratio of specific heats 
was fixed at 1.2, and the constants for Sutherland’s law were set ac- 
cording to values appropriate for steam. These values were chosen 
as representative of the combusting gaseous hydrogen and gaseous 
oxygen flow with an oxidizer-to-fuel ratio of six. The values were 
obtained by examining a series of one-dimensional equilibrium 
calculations 22 for the chamber conditions. The total temperature in 
the chamber was observed to increase slightly with increasing total 
pressure. The value of the gas constant was also observed to change 
a few percent over the pressure range of interest. This translates into 
a slightly different value of reference Reynolds number for each of 
the three chamber pressures considered. Thus the chamber condi- 
tions were specified as shown in Table 2. The Prandtl number was 
fixed at 0.6. The compositions of the primary and secondary flows 
were assumed to be equivalent. This was a constraint set by utilizing 
a single-species, calorically perfect-gas code, but nevertheless befits 
the assumption of a gas-generator-cycle engine. The secondary flow 
was injected at a constant total temperature of 1600°R. 

D. Mesh Generation 

Mesh generation was executed with Gridgen vl t software. 23 Sev- 
eral rules were observed for all of the different mesh geometries to 
ensure that similar spatial resolution was achieved for all cases: 

1 ) Define the grid point clustering toward all no-slip surfaces with 
a near-wall spacing of 3.0 x 10 4 (/ ref ). 

2) Define maximum spacing away from surface not to exceed 
0.25 (/ ref ). 

3) Define maximum stretching ratio of 1.2 between adjacent cells 
in any one dimension. 

4) Utilize a three-block mesh topology with point-to-point in- 
terblock boundaries. 

The price of adequate resolution was a large number of mesh 
points for the largest ducts (90,701). The important issue of doc- 
umenting the mesh independence is addressed in the following 
section. 

V. Statistical DOE 

A. DOE and RSM 

DOE is grounded in statistics and, thus, takes a radical departure 
from the typical Edisonian approach to experimentation. In the field 


of DOE, the primary objective is to effectively model the response 
of an experiment, as a function of the independent variables. Once 
the experiments are complete, an empirical model is constructed and 
can be optimized to meet overall objectives of reliability and perfor- 
mance. One advantage of this approach lies in the fact that optimum 
performance is not limited to the finite set of experimental responses 
gathered during experimentation. Rather, an optimum response usu- 
ally exists in the design space at a region between the completed ex- 
periments. Prior to optimization, the final step in model construction 
involves validation of the empirical model. Successful prediction of 
a few additional experiments will build confidence in the empirical 
model. Of course, uncertainty analysis governs the reliability as- 
sessment of the empirical model. Therefore, reproducibility of the 
response must be quantified. 

DOE strategies also have the advantage of minimizing the num- 
ber of experiments, here CFD simulations, based on the complexity 
of the anticipated relationships between the response and the inde- 
pendent variables under investigation. In our specific situation the 
design complied with several other a priori restrictions. 

1 ) The initial subset of nine CFD solutions should result in the 
formation of a multiple linear regression (MLR) model (hyperplane) 
for the purposes of comparing turbulence models and screening the 
independent variables to ensure their ranges were properly chosen. 

2) The entire set of CFD solutions should result in the formation 
of an RSM model with linear, curvilinear, and two-way interaction 
terms. 

3) The total number of required CFD solutions should result in a 
practical number of parametric experiments (i.e., <40). 

The final design selected was a 36-case D-optimal 24 RSM model 
with an embedded multiple linear regression model. In other words 
the first nine CFD solutions permitted the estimation of a simple sta- 
tistical model for general predictions of nozzle performance, and the 
full 36 solutions supported a higher-order model for accurate pre- 
dictions of nozzle performance. The higher-order model contained 
an intercept [c 0 ], six linear terms [cjjti, 02 * 2 , ■ • • . c 6 .r 6 ], 15 bilinear 
or interaction terms [c^x^xi, c^x\x^ t . . . , 021 at 5 A 6 ], and six curvilin- 
ear terms [c 22 *i , C 23 JC 2 , . . . , cnx\\. In practice, only the coefficients 
(c 0 , c 1 , . . . , C 27 ), which are statistically significant (statistically dif- 
ferent from zero), are retained in the final solution. Thus, DOE re- 
duced the full parametric study (or full factorial) of 729 cases down 
to just 36. The final RSM model, once checked for adequacy of fit to 
the response data, was capable of predicting nozzle performance at 
any point within our six-dimensional design space without having 
to execute additional CFD simulations. The D-optimal design con- 
struction, MLR analyses, and RSM model analyses were executed 
with RS/l software . 25 Detailed specification of the six independent 
variables for each case is presented in Table 3. Note that case 9 was 
the design centerpoint; each variable was set at the midlevel value. 

Great care was taken to reduce potential sources of computa- 
tional bias. Grid dependence, incomplete convergence, and numer- 
ical instability are all possible sources of bias in the response data 
that could potentially lead to false conclusions. Furthermore, unlike 
physical experimentation, repeating a CFD experiment should be 
exactly reproducible within machine accuracy. Hence, the typical 
MLR strategy of fitting a statistical model down to the level of ex- 
perimental noise had to be abandoned. Instead, we fit the statistical 
model down to a level of uncertainty no less than ± 1 %, to avoid 
overfitting the data. In other words, we expected to calculate an ac- 
curate representation of the governing model equations and report a 
response to within approximately ±1% uncertainty. In the next few 
sections we document the approach used to control these sources of 
computational bias and defend this assumed uncertainty estimate. 

B. DOE and Mesh Dependence 

A detailed mesh-convergence analysis was performed for the de- 
sign centerpoint, case 9. This case was intentionally included as one 
of the 36 cases tested (from a possible matrix of 3 6 possibilities). It 
was a good candidate for detailed analysis because of its symmet- 
ric placement in our six-dimensional design space. The standard 
discretization of case 9 resulted in 56,952 points in the computa- 
tional domain. A doubling of the grid in each direction results in * 
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Table 3 Thirty-six case experimental matrix: 
independent variable specification 


Case 


Independent variable settings 



(P c , psi) 

xi 

(, m s x 100 %/riip) 

-*3 

{Ay/A*) 

X 4 

(L/Dy) 

*5 

{Ad Ay) 

X6 

(A r /A*) 

1 

300 

8 

39 

2 

2 

4 

2 

300 

0 

39 

5 

2 

20 

3 

1200 

0 

200 

2 

2 

4 

4 

300 

0 

200 

5 

1 

4 

5 

300 

8 

200 

2 

1 

20 

6 

1200 

8 

200 

5 

2 

20 

7 

1200 

8 

39 

5 

1 

4 

8 

1200 

0 

39 

2 

1 

20 

9 

750 

4 

120 

3.5 

1.5 

12 

10 

300 

0 

200 

2 

2 

20 

11 

1200 

8 

39 

2 

2 

20 

12 

1200 

0 

200 

5 

1 

4 

13 

1200 

0 

39 

5 

2 

20 

14 

300 

0 

200 

5 

2 

4 

15 

1200 

8 

200 

5 

1 

20 

16 

300 

8 

200 

2 

2 

4 

17 

300 

8 

39 

5 

2 

20 

18 

1200 

U 

39 

2 

1 

4 

19 

300 

0 

200 

5 

1 

20 

20 

300 

8 

200 

5 

1 

4 

21 

1200 

8 

39 

5 

2 

4 

22 

1200 

0 

200 

2 

1 

20 

23 

750 

0 

39 

5 

1 

4 

24 

750 

4 

200 

2 

1 

4 

25 

300 

4 

39 

5 

1.5 

4 

26 

750 

8 

120 

2 

1.5 

4 

27 

1200 

0 

200 

5 

1.5 

12 

28 

300 

4 

120 

5 

1 

12 

29 

300 

8 

39 

3.5 

1 

4 

30 

300 

0 

120 

2 

1 

4 

31 

300 

4 

39 

2 

1 

20 

32 

1200 

8 

200 

2 

1 

12 

33 

300 

0 

39 

3.5 

1.5 

20 

34 

750 

0 

39 

2 

2 

12 

35 

1200 

8 

39 

5 

l 

20 

36 

1200 

4 

120 

3.5 

2 

4 


Table 4 Mesh dependence of 7 sp for case 9 


Parameter 

Standard mesh 

Fine mesh 

# mesh point 

56,952 

226,281 

/sp, s 

418.4 

4169 

A% 

0.346% 

GCI (p — 1.5) 

1.6% 


quadrupling of the total points, to 226,281. The integral quantities 
changed very little, as evidenced by the data of Table 4. 

The mcsh-dependence data were presented in terms of a net differ- 
ence and the standard grid-convergence index (GC1). 26 The unique 
advantage of reporting the GCI value lies in the uniform manner of 
presenting dependent variable changes. The GCI is not a bound on 
the error; rather it is . .a reasonable error band, in the flavor of a 
statistician’s 2<r range or an experimentalist’s 20: 1 odds. . 25 The 

GCI is calculated with the following formula: 


GCI = 


3(A%)<*^ 
(at* - 1 ) 


(4) 


where a is the grid-refinement ratio and p is the numerical order of 
accuracy of the CFD solver. For our study, a = 2 and conservatively 
estimating the spatial accuracy of the CFD between first and second 
order, we set P — 1.5. The GCI value was 1.6%. Thus, we assume 
that the standard approach to mesh generation just outlined has 
removed any significant mesh dependence of the / sp results upon 
which our statistical model is based. 

The local flow properties of wall static pressure and centerline 
Mach number, shown in Fig. 2, revealed noticeable mesh depen- 


226,952 pt. mesh 

56,952 pt. mesh 

Pw»ij 226,281 pt. mesh 

p W4lJ 56,952 pt. mesh 



Axial Distance From Throat (in) 


Fig. 2 Mesh dependence: centerline Mach number and mixer-ejector 
surface pressure for case 9. 


3 

S 



Axial Distance From Throat (in) 


Fig. 3 Mass-flux ratio vs axial location, plotted from the rocket exit 
plane to the mixer-ejector exit plane for case 9. 


dence in the solution of case 9. This does not, however, contradict 
the data of Table 4. One would expect integral properties to exhibit 
less sensitivity to mesh refinement than local properties, and thus 
reach an asymptotic value earlier in a mesh-refinement study. This 
is particularly true for local properties in the neighborhood of strong 
gradients such as shock and contact surfaces. We have included the 
data of Fig. 2 to highlight an important observation: different aspects 
of a given CFD simulation reach mesh independence at different lev- 
els of mesh refinement. Mesh-refinement analysis must be carried 
out upon the specific response data of interest. Clearly a response 
surface model build upon local centerline Mach-number data would 
require further mesh refinement to avoid an unacceptable level of 
bias/uncertainty caused by inadequate spatial resolution. 


C. DOE and Iterative Convergence 

At the outset, we addressed the issue of the reported values of 
/ sp and their sensitivity to incomplete iterative convergence. Often 
overlooked, this issue proved to be nontrivial when considering DOE 
and CFD. The residuals of mass flow and thrust converged at least 
four orders of magnitude. Furthermore, we demanded that the value 
of / sp should not change more than 0.15% over the last 10,000 it- 
erations at a Courant number of 0.5. This criterion required several 
of the longest domains to be run in excess of 50,000 iterations, de- 
pending upon the initial conditions prescribed. The extremely long 
iteration counts were necessary to resolve the widely varying wave 
speeds present in a domain with large regions of low-speed flow 
present alongside hypersonic jet flow. Figure 3 illustrated the small 
variation of integrated mass flux as a function of axial location for 
the converged result of design centerpoint, case 9. The mass-flux 
ratio of Fig. 3 was simply the local integrated mass flux, normalized 
by the average integrated mass flux. We assume that the preceding 
approach for convergence analysis has removed any significant iter- 
ative dependence of the / sp results upon which our statistical model 
is based. 


D. DOE and Artificial Dissipation 

For most of the 36 cases, the coefficients of second- and fourth- 
order artificial viscosity were fixed at £ and respectively. 
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However, for some of the configurations involving the largest pri- 
mary rocket nozzles, numerical stability became an issue. Conver- 
gence of the exit mass flux and thrust were oscillating such that / sr 
varied as much as ± 1 .5% instead of converging toward machine ac- 
curacy. We found that specifying the coefficient of fourth-order dif- 
fusion at a value of ~ was sufficient to stabilize the flowfield without 
altering the / sp integrals significantly. All of the cases run with the el- 
evated fourth-order artificial diffusion had / sp values that converged 
within the band of oscillation already observed. No adjustment to 
the second-order artificial diffusion coefficient was needed for any 
of the 36 cases. 

To be complete about assessing the effects of elevated fourth- 
order artificial diffusion, we decided to examine a case that did not 
suffer from numerical instability. The design centerpoint was a good 
candidate for the reasons just mentioned. When the CFD simulation 
of case 9 was repeated with the higher coefficient of fourth-order 
artificial diffusion, the / sp value changed from a converged value of 
418.36 to 419.65 or 0.30%. Thus, we assume that the modification 
of fourth-order diffusion just outlined has removed the uncertainty 
associated with numerical stability without compromising the / sp 
results upon which our statistical model is based. 

VI. Results 

A. RBCC Nozzle Flowfield 

A few dominant flow features were common to all RBCC noz- 
zles we have examined. A brief examination of these phenomena 
provided for a more complete interpretation of the statistics. Figure 
4 is a composite picture of the design centerpoint, case 9. The upper 
half of the figure is a numerical schlieren image, and the lower half 
is a plot of iso-Mach contours. The numerical schlieren image was 
based upon a transformation of the calculated density gradient. 27 It 
was particularly helpful in revealing the presence of expansion fans, 
shock waves, and shear layers. 

In the following discussion, we will refer to the expansion within 
the RBCC nozzle in three stages, shown in Fig. 4. The first or primary 
stage of expansion occurred within the contoured rocket nozzle. 
The secon I or free expansion, occurred between the rocket exit and 
mixer-ej 'ctor wall. The third, or duct expansion, occurred within 
the diverg ig mixer-ejector. 

Two dif erent compression phenomena were observed. A con- 
tinuous ser es of weak compression waves were visible within all 
three of the Rao-optimized 28 rocket nozzles studied. Also note that 
a plume shock formed downstream of the rocket nozzle exit and 
always reflected off of the duct wall. For elongated mixer-ejector 
geometries the plume shock was seen to reflect several times be- 
tween the duct wall and centerline. 

B. Flow Visualization 

The primary expansion was governed by the rocket nozzle ge- 
ometry. This nozzle contour was optimized with the Rao algorithm. 


Plume Coalescing 

Shock Compression Weak Nozzle 



M 



03 U 2.7 34 44 53 6A 


A perfect nozzle contour balanced the expansion and compression 
waves generated at the cost of excessive nozzle length. A Rao con- 
tour shortened the nozzle length while maximizing the thrust gen- 
erated. The cost of this optimization was the generation of weak 
compression waves from the contoured nozzle. These compression 
waves were present within all three of our Rao-contoured nozzles 
and were visible in Fig. 4. The compression waves originated at the 
nozzle wall and coalesced at some point downstream of the rocket 
nozzle exit plane. It appeared that these weak waves had begun 
to coalesce along the upstream edge of the compression and fully 
coalesced after reflecting off of the centerline. 

The free expansion revealed a structure similar to an underex- 
panded supersonic free jet. A high-speed shear layer surrounded the 
primary plume shock. The shape of this high-speed shear layer was 
a function of the rocket nozzle exit pressure and the static pressure 
of the annular cavity region. This high-speed shear layer defined the 
boundary between the fast moving plume flow and the slow (or stag- 
nant) cavity flow. Because the shear layer was a constant pressure 
process, its location was a function of the cavity pressure. 

The static pressure of the annular cavity was dependent upon the 
cross-sectional area and the magnitude of secondary flow injected. 
Secondary flow injection also had an important effect in the stagna- 
tion region of the shear layer. Consider the pictures of cases 20 and 4 
shown in Fig. 5. The only difference between these two simulations 
was the presence of 8% secondary flow injected into the annular 
cavity. Note the different locations of the primary shear layers and 
shock waves. An examination of the peak static pressure along the 
duct wall (Fig. 6) revealed that the primary shock wave impinged 
further downstream in the presence of secondary flow injection. The 
impingement point of case 20 was nearly twice that of case 4, when 
measured from the rocket exit plane. This relates directly to the 
compression necessary to tum the flow toward the axial direction. 

Figure 5 also dramatically illustrated the impact that injected sec- 
ondary flow had upon the development of streamline patterns inside 



Fig. 5 Particle traces: gray represents the secondary flow streaklines; 
black represents the primary flow streaklines for a) no injected flow 
(case 4) and b) 8% injected flow (case 20). 



Axial Distance from Rocket Exit (in) 


Fig. 4 Composite image of design centerpoint, case 9: a) numerical 
schlieren image and b) computed Mach-number contours. 


Fig. 6 Mixer-ejector surface pressure comparison; notice the effect of 
secondary flow upon the peak pressure (shock impingement). 
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Table 5 Tabulation of the / sp data from all of the CFD simulations 


Case 

/ MAX 
sp 

Calculated / s c p ro 

Calculated 

r,f° 
-LLsl 

k-E 

SA 

k~£y % 

SA, % 

] 

451.57 

410.42 

399.54 

90.89 

88.48 

2 

451.57 

427.06 

424.33 

94.57 

93.97 

3 

480.04 

415.88 

415.94 

86.63 

86.65 

4 

465.06 

373.85 

361.73 

80.39 

77.78 

5 

46506 

425.22 

411.29 

91.43 

88.44 

6 

480.04 

437.21 

430.64 

91.08 

89.71 

7 

445.44 

375.30 

374.07 

84.25 

83.98 

8 

445.44 

418.70 

416.89 

94.00 

93.59 

9 

468.03 

423.38 

418.36 

90.46 

89.39 

10 

473.31 

— 

433.80 

— 

91.65 

11 

457.99 

— 

425.32 

— 

92.87 

12 

471.67 

— 

375.90 

— 

79.69 

13 

457.99 

— 

432.92 

— 

94.53 

14 

473.31 

— 

393.85 

— 

83.21 

15 

471.67 

— . — 

415.28 

— 

88.04 

16 

473.31 

— 

407.86 

— 

86.17 

17 

451.57 

— 

412.17 

— 

91.28 

18 

445.44 

— 

376.14 

— 

84.44 

19 

465.06 

— 

403.11 

— 

86.68 

20 

465.06 

— 

372.02 

— 

79.99 

21 

457.99 

— 

399.24 

— 

87.17 

22 

471.67 

— 

419.26 

— 

88.89 

23 

443.34 

— 

369.83 

— 

83.42 

24 

469.45 

— 

381.15 

— 

81.19 

25 

446.71 

— 

384.22 

— 

86.01 

26 

468.03 

— 

404.67 

— 

86.46 

27 

476.73 

— 

415.05 

— 

87.06 

28 

458.00 

— 

397.36 

— 

86.76 

29 

439.19 

— 

377.51 

— 

85.96 

30 

458.00 

— 

369.02 

— 

80.57 

31 

439.19 

— 

399.86 

— 

91.04 

32 

471.67 

— 

415.49 

— 

88.09 

33 

446.71 

— 

419.92 

— 

94.00 

34 

455.83 

— 

427.36 

— 

93.76 

35 

445.44 

— 

401.40 

— 

90.11 

36 

473.98 

— 

404.55 

— 

85.35 


the annular cavity region. If the cavity was closed (m s /m p = 0), 
then the secondary flowfield was essentially a supersonic driven 
cavity. Case 4 developed two large recirculating zones and several 
smaller ones. A substantial shear layer between the two large ed- 
dies was visible in many of the schlieren images of closed cavities. 
Case 20, however, developed one small eddy near the primary shock 
impingement. Otherwise, the secondary flow was sent downstream 
along the duct wall. 

The third stage of expansion occurred inside the mixer-ejector 
only when divergence was present (A 6 / A3 > 0). Oblique shock 
waves often reflected at the jet centerline and the mixer-ejector 
wall. The shock- wave/boundary -layer interaction often resulted in 
flow separation and reattachment. For these cases an oblique shock 
formed off the front and back sides of the larger separation bub- 
bles. This complex behavior is well known to occur for hypersonic 
shock- wave/boundary-layer interactions. 29 

C. Calculated Nozzle Efficiency 

We shall begin by presenting the tabulation of the data, from all 

of the CFD simulations, in Table 5 . The first nine cases were executed 
twice, once with each of the turbulence models under consideration. 
This allowed an independent regression analysis for each turbulence 
model. Analysis of variance techniques were applied to determine if 
the differences between the two turbulence models were statistically 
significant. Cases 10-36 were calculated with the S-A turbulence 
model only. The entire 36 case data set using S-A was used to build 
the full RSM model described next. 

D. Multiple Linear Regression Analysis: Cases 1-9 

Both MLR analyses of the two turbulence models agreed upon 
the significance and magnitude of four linear effects. Primary rocket 
area ratio and mixer-ejector exit area exhibited a positive linear 


relationship with Mixer-ejector inlet area and mixer-ejector 
length/diameter ratio exhibited a negative linear relationship with 
rj lsp . However, the secondary flow and chamber pressure did not 
appear to have a significant effect upon nozzle efficiency for either 
turbulence model. 

The analysis of variance revealed a small, albeit statistically, sig- 
nificant difference between the two turbulence models’ predictions. 
This difference was in the magnitude of the intercept value; k-£ 
model coefficient was approximately 1% greater. This simply indi- 
cated that the k-e model tended to predict a slightly higher level 
of nozzle efficiency. We interpreted this to indicate that the choice 
between these two turbulence models would not affect the relative 
prediction of the regression, but would slightly influence the ab- 
solute prediction. This can be confirmed by briefly examining the 
slightly higher efficiencies shown for k-e in Table 5. A lengthy 
description of these linear regression results has been omitted in 
favor of concentrating upon the more complete RSM results. The 
interested reader is referred to a companion paper. 18 

The preceding MLR analysis was available to us before the com- 
pletion of cases 10-36. As a result, we chose to use the S-A model 
instead of the k-e model to complete the study for the following rea- 
sons: 1) faster iterative convergence, 2) more robust performance for 
our configurations, and 3) statistically similar efficiency predictions. 
Therefore, all solutions discussed from here on will refer to the S-A 
model results. 

At this point the MLR analysis produced a model, which revealed 
that nozzle efficiency varied linearly with four different paramet- 
ric effects. However, the MLR analysis had to assume that neither 
curvature nor interaction existed in this regression. This was an in- 
trinsic assumption from the chosen nine-expenment DOE. We felt 
that these effects might actually be present because so much of fluid 
analysis is truly nonlinear. One hint of the nonlinearity lurking in the 
results was given by comparing the MLR-predicted efficiencies with 
the actual CFD values for the first nine cases. This revealed that the 
linear model fit the CFD data at the edges of design space better than 
at the center. This can imply curvature in the data, not accounted for 
in the simple regression model. Another concern involved the fact 
that we had fit a five-term model (intercept and four linear terms) 
using nine CFD simulations. This leaves us with just four extra de- 
grees of freedom in the regression error, a statistical entity used to 
access model adequacy. By completing the full 36 cases, we directly 
addressed the question of nonlinearity while gaining more degrees 
of freedom in the statistical model’s error term. 

E. Response Surface Model Analysis: Cases 1-36 

As expected, the response surface model revealed many more sig- 
nificant effects than the simpler linear regression analysis. In total, 
13 separate model coefficients were considered significant in the 
RSM analysis [Statistically, the model coefficients are either sig- 
nificantly different than zero (i.e., greater than 95% significant), or 
assumed to be zero.]; the model coefficient values and significance 
levels are given in Table 6. For the purposes of the following dis- 
cussion, the model terms of the form Jt, will be referred to as linear, 
XjXj as bilinear, and xf as curvilinear. The statistical model to pre- 
dict RBCC system efficiency was defined by Table 1, Eq. (1). and 
Table 6. It is given in Eq. (5): 

^rsm = (8 9325) + (0.3503)ii - (1.6150)i 3 - (0.6976)* 4 
+ (L6077)*5 + (3.2062)Jc 6 -F (0.3344 )Jc,Jc 3 
+ (0.3626)jc 2 i 3 “ (0.4209)i 2 ^5 - (0.7228)^^6 

+ (0.2830).ij* 5 - (0.406 1 )V„ - (1 -5860)jc| (5) 

An approximate prediction error and specified confidence inter- 
val can be calculated for this RSM. This error estimate is based 
upon the root-mean-squared error of the regression and Student’s 
t-distribution. The 95% confidence predictive uncertainty associated 
with the RSM model was approximately (±1.1%). 

Although Eq. (5) provided a means for quantitative examination 
of our RSM model, graphical interpretation was far more useful for 
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Table 6 RSM model coefficient values 
and significance 


RSM model 
term 

Coefficient 

Significance, % 

1 

89.325 

99.99 

*1 

0.3503 

99.16 

* 2 

0.0 

— 

*3 

-1.6150 

99 99 

*4 

-0.6976 

99.99 

*5 

1.6077 

99 99 

*6 

3.2062 

99.99 

*1*2 

0.0 

— 

*1*3 

0.3344 

98.46 

*1*4 

0.0 

— 

*1*5 

0.0 

— 

*1*6 

0.0 

— 

*2*3 

0.3626 

98.89 

*2*4 

0.0 

— 

*2*5 

-0.4209 

99.53 

*2*6 

-0.7228 

99.99 

*3*4 

0.0 

— 

*3*5 

0.2830 

95.53 

*3*6 

0.0 

— 

*4*S 

0.0 

— 

*4*6 

0.0 

— 

*5*6 

-0.4061 

99.45 

*? 

0.0 

— 

*? 

0.0 

— 

*? 

0.0 

— 

*4 

0.0 

— 

£ 2 
*5 

0.0 

— 

A 

*6 

-1.5860 

99.99 



Fig. 7 Efficiency ry as a function of mixer-ejector inlet area ratio 
(A 3 /A*) and rocket area ratio (A r /A*). 

exploring the results. Each of the following sections begins with a 
three-dimensional carpet plot of the efficiency, followed by a dis- 
cussion of the physical implications. For each carpet plot, efficiency 
was portrayed along the vertical axis as a function of two indepen- 
dent variables. Unless otherwise noted, the remaining independent 
variables are set to their respective midlevel values. 

I. Annular Cavity Size 

The carpet plot of efficiency ry as a function of rocket area ra- 

tio (A r /A*) and mixer-ejector inlet area ratio (A 3 /A*) revealed 

a large part of our story. This carpet plot, shown in Fig. 7, por- 

trayed efficiency as a function of two geometric variables. Thus 

each point on the surface corresponded to a slightly different ge- 
ometry, with the limiting configurations shown in two-dimensional 

silhouette at the four comers. This surface was curvilinear because 


of the quadratic influence of rocket area ratio (i.e., term x\). Also 
note that the four comer points are coplanar, evidence that no inter- 
action occurs between these variables. 

Examination of Fig. 7 highlighted two important features: 1) an 
increase in rocket area ratio corresponded to an increase in perfor- 
mance, and 2) a decrease in mixer-ejector inlet area ratio led to an 
increase in performance. Consequently, the maximum performance 
shown on this surface occurred for the configuration (A r /A* =20) 
and (A)/ A* = 50). The physical implication of these two observa- 
tions was consistent: reduce the size of the annular cavity region to 
improve the efficiency. Recall that reduction of the annular cavity 
corresponds to a reduction of the free-expansion portion of the flow. 
These observations lead to the suspicion that free expansion across 
the annular cavity was an inefficient process, when compared to the 
attached expansions within the Rao nozzle and the mixer-ejector. 
Further analysis of the CFD data confirmed this suspicion. 

To measure the thrust achieved during the free expansion, we in- 
tegrated the static pressure in the annular cavity across the upstream 
wail. However, this method of calculating the free-expansion thrust 
was strictly valid for the subset of cases with straight mixer-ejectors 
and no injected secondary flow. The calculated free-expansion thrust 
was then compared with the isentropic value expected from an area 
increase of (A r /A*) to (A 3 /A*). The ratio of these two values of 
thrust can be expressed as the efficiency of free expansion. Topical 
free-expansion efficiency values were between 20-30%, depending 
on the specific geometry. This was in sharp contrast to the Rao- 
optimized rocket nozzle expansions, which registered above 95% 
efficient for all cases. Static pressure observed within the cavity was 
essentially constant behind the high-speed shear layer. This obser- 
vation agreed with the data of Fortini. I. * * * * * 7 

We believe that the losses associated with the large free expansion 
are caused primarily by overexpansion. We have already noted that 
the cavity pressure drops below the value predicted by isentropic 
flow, which caused the expanding flow to exit the rocket nozzle at 
a steeper angle than anticipated. This required a stronger oblique 
shock structure to turn the flow toward the axial direction, and ulti- 
mately led to a more intense stagnation flow at the duct wall. Exam- 
ine the flowfield of case 4 depicted in Fig. 5. Notice how steeply the 
plume expanded in the radial direction, downstream of the rocket 
nozzle exit. 

2. Mixer-Ejector Divergence 

The carpet plot of efficiency t] as a function of mixer-ejector 
inlet area ratio (A 3 /A*) and exit area ratio (A 6 /A 3 ) revealed an- 
other important part of our study. Figure 8 clearly shows that in- 
creasing mixer-ejector exit area ratio substantially improved the 



Fig. 8 Efficiency r? as a function of mixer-ejector inlet area ratio 
(A 3 /A*) and mixer-ejector exit area ratio (A^A)). 
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Case #15 



Numerical schlieren image 
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Mixer-ejector surface- pressure comparisons 

Fig. 9 Effect of mixer-ejector divergence. 

performance. Again, we see that decreasing the mixer-ejector inlet 
area improved performance. For the configurations shown in Fig. 8, 
the maximum efficiency occurred in the comer where (A 3 /A* = 50) 
and (A 6 /A 3 = 2). 

Increasing (A 6 /A 3 ) to a value greater than unity was tantamount 
to introducing divergence into the mixer-ejector geometry. This 
divergence improved the efficiency for two reasons. Divergence 
caused an efficient attached-expansion to occur within the mixer- 
ejector. Divergence also caused the plume shock to reflect at a shal- 
lower angle, which delayed the onset of boundary-layer separa- 
tion downstream. These effects were visible in the comparison of 
cases 6 and 15 in Fig. 9. The physical implications were consis- 
tent with that just discussed. Namely, the efficiency was maximized 
when the influence of free-expansion losses was reduced as much as 
possible. 

The slope of the response surface of Fig. 8, with respect to 
(A 6 /A 3 ), was a function of (A 3 /A*). In other words, the slope 
of this carpet plot [9rj(jc 3 , Jc 5 ) / 3 Jc 5 ]jf t = o depended upon the 
value of (x 3 ). This was caused by the interaction between these 
two variables, or term Jt 3 x 5 of Eq. (5). In fact, the slope of the six- 
dimensional response surface with respect to x 5 depended upon the 
values of Jc 2 , i 3 , jc 6 , caused by three separate bilinear interactions. 
Despite the seeming complexity of variable interactions, the follow- 
ing observations can be made: 1) the surface gradients of Fig. 8 were 
largest near the configuration with minimum efficiency, and 2) the 
free-expansion losses were maximized for the same configuration 
because of the large mixer-ejector inlet ratio and no mixer-ejector 
divergence. Thus, the variable interaction corroborated our earlier 
evidence. 

3. Rocket Area Ratio and Mixer-Ejector Divergence 

The carpet plot of efficiency rj as a function of rocket area ratio 
(A r /A*) and mixer-ejector exit area ratio (A 6 /A 3 ) revealed the third 
major part of the results. Figure 10 confirmed what we have already 
observed from Figs. 7 and 8: 1) increasing (A r /A*) improved effi- 
ciency and 2) increasing (A6/A 3 ) also improved efficiency. Again, 
the quadratic effect of rocket area ratio upon performance was re- 
vealed in the curvilinear surface plot. However, the two quadratic 
curves that bound this surface were different, because of the variable 
interaction. Thus, this surface plot has curvature because of the term 
x\ and twist caused by the bilinear term x 5 x 6 . 

We have already discussed the influence of both (A r /A*) and 
(A 6 /A 3 ) upon the efficiency. However, the variable interaction .£ 5 X 6 
was not visible from Figs. 7 and 8. We can see from Fig. 10 that the 
configuration with the largest annular cavity and no mixer-ejector 
divergence offered the least efficient performance. Consequently, 
this configuration stood to gain the most by increasing ( A r /A *) 
and/or increasing (A 6 /A 3 ). Thus, the largest surface gradients are 
located in this comer of the carpet plot and gradually decrease in 
magnitude until the most efficient configuration was achieved. The 
variable interaction x 5 x 6 was evidence of this phenomenon. 



Fig. 10 Efficiency tj as a function of mixer-ejector exit area ratio 
(A$M 3 ) and rocket area ratio (, A r /A *). 



Fig. 11 Efficiency rj as a function of mixer-ejector length/diameter 
ratio (ZVJ) 3 ) and mixer-ejector exit area ratio (A 6 /A 3 ). 


4. Mixer- Ejector Length 

The carpet plot of efficiency as a function of mixer-ejector exit 
area (A 6 /A 3 ) and length/inlet diameter (L/D 3 ) ratios was quite 
interesting. Figure 1 1 clearly shows the small influence that mixer- 
ejector length had upon efficiency. Extending (L/D 3 ) from two 
to five resulted in a few percent reduction of efficiency. We at- 
tributed this to the drag resulting from increased wetted surface 
area within the mixer-ejector. Notice that for a given (Z,/D 3 ), in- 
creasing (A 6 /A 3 ) resulted in a net gain in efficiency, despite the 
increase in wetted area. The key point to emphasize here is that by 
increasing the duct exit area we diluted the free-expansion losses by 
including an additional attached expansion within the mixer-ejector. 
On the other hand, a simple extension of the mixer-ejector geometry 
merely added drag to the system, without any increase in thrust. The 
response surface of Fig. 1 1 was a simple plane without any variable 
interaction. 

5. Secondary Flow Injection 

Recall Fig. 7, which displayed the efficiency as a function of 
mixer-ejector inlet area ratio (A 3 / A*) and rocket area ratio ( A r /A *). 
We observed that maximizing (A r /A*) was important along with 
minimizing (A 3 /A*). Secondary flow (mjm p ) had an interacting 
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□ m s / rh p = 0% 



Fig. 12 Efficiency rj as a function of mixer-ejector inlet area ratio 
(Aj/A*) and rocket area ratio (A r /A*): interaction of secondary flow 
(mjm p ). 

effect with both of these parameters, as seen from the terms x 2 x 3 
and x 2 X(, in E<T (5). 

Figure 12 displays the effect of three different values of (m s /m p ). 
The gray surface, similar to that plotted in Fig. 7, corresponded 
to the midlevel of (m y /m p ~ 4%); the white surface corre- 
sponded to ( m s /m p =0%); and the dark surface corresponded to 
(m s /m p = 8%). Plotting all three surfaces together highlighted the 
influence of secondary flow upon the effects of these two variables. 
The extrema of performance were altered by secondary flow injec- 
tion, but in opposite ways. The performance of the most efficient 
configurafion was maximized without injecting any secondary flow. 
However, : icreasing the secondary flow injection to 8% maximized 
the perfori ance of the least efficient configuration. Secondary flow 
injection ii to the annular cavity must have helped to relieve the 
free-expan* ion losses, but at an overall system cost. 

We suspect that the benefit of secondary flow derived from the 
increased static pressure within the cavity. This had the effect of 
reducing the free expansion, thereby reducing both the expansion 
angle and the strength of the reflected plume shock. There was also 
a direct contribution to stream thrust from secondary flow injec- 
tion, although the / sp of this secondary flow was dependent upon 
the precise configuration. We suspect that injection of secondary 
flow, with a lower total pressure, into the mixer-ejector was a com- 
peting loss mechanism. The cost-benefit analysis of secondary flow 
injection appeared to rest on the significance of the free-expansion 
losses, which were geometry specific. We also suspect that a dif- 
ferent composition of the secondary flow may alter the balance of 
these competing effects. 

Secondary flow had an interactive effect with mixer-ejector exit 
area ratio also. This effect was quite similar to that shown in 
Fig. 12, in that the extrema of performance were affected. From 
Figs. 8 and 10, we have shown that performance increased when the 
mixer-ejector exit area was maximized. We surmised that the free- 
expansion losses were most significant for mixer-ejectors without 
divergence (A6/A 3 = 1). Thus we might expect that maximizing 
secondary flow injection would improve performance for mixer- 
ejectors with minimum divergence and likewise reduce performance 
for mixer-ejectors with maximum divergence. This was precisely 
what occurred. 

F. Results Summary 

1 ) Attached expansions within the rocket nozzle and mixer-ejector 
duct were far more efficient than free-expansion across the annular 
cavity. 


Table 7 Confirmation runs to examine the bilinear interaction 
(£14*2*5) (Note that xu * 3 » * 4 , and x 6 are fixed at their mean values.) 


Case 

*2 (ms/nip), % 

*5 (At/ As) 

rjCPD % 

'* p 

RSM % 
'*P 

A, % 

37 

0 

2 

89.94 

91.35 

-1.41 

38 

0 

1 

86.98 

87.30 

-0.32 

39 

8 

2 

89.65 

90 51 

-0.86 

40 

8 

1 

87.58 

88.14 

-0.56 

Table 8 Confirmation runs to examine the bilinear interaction 
(ci 5 * 2 * 6 > (Note that jcj, x 3 ,x 4 , and x$ are fixed at their mean values.) 

Case 

*2 ("ly/mp), % 

*6 {A r t A*) 

i£ FD ,% 

/sp 

„? sm ,% 

/sp 

A, % 

41 

0 

20 

91.29 

91.67 

-0.38 

42 

0 

4 

83.46 

83.81 

-0.35 

43 

8 

20 

89.69 

90.22 

-0.53 

44 

8 

4 

85.39 

85.26 

0.13 


2) Nozzle efficiency was improved when the size of the annular 
cavity was reduced by increasing the rocket exit area ratio (A r / A*) 
and reducing the mixer-ejector inlet area (A 3 /A*). 

3) Mixer-ejector expansion, present when (A 6 /A 3 > 1), generally 
improved efficiency by efficiently expanding the flow within the 
mixer-ejector and reducing the losses associated with the reflecting 
plume shock. 

4) Nozzle efficiency was inversely proportional to the overall 
length of the mixer-ejector duct; this was a relatively minor effect. 

5) Secondary flow injection (m s /m p ) could mitigate free- 
expansion losses, although at some overall cost, through variable 
interactions: relatively inefficient configurations, with large free- 
expansion regions, benefited from maximizing (m s /m p ) and rela- 
tively efficient configurations, with small free-expansion regions, 
benefited from minimizing ( m s /m p ). 

6) Statistically significant bilinear and curvature terms under- 
scored the importance of executing the full 36-case D-optimal de- 
sign versus the simple nine-case MLR design. 

VII. RSM Confirmation 

The final step in DOE involved the execution of confirmation 
tests. These additional computations tested the overall ability of the 
RSM model to predict RBCC nozzle efficiency integrated from CFD 
simulation. However, the eight new cases were chosen to specifi- 
cally assess the validity of the two largest bilinear interactions x 2 x 5 
and x 2 x 6 . Cases 37-40 were configured to address the former bi- 
linear interaction, and cases 41-44 addressed the latter, as shown 
in Tables 7 and 8. All but case 37 were well within the approxi- 
mate uncertainty band of (±1%). Case 37 results lay just outside 
the uncertainty band. 

VIII. Conclusions 

The present computational effort represents an essential first step 
in assessing the performance of RBCC systems during all-rocket 
mode operation. We have utilized a statistical design of experiments 
approach to quantify the influence of four geometric and two flow- 
field variables using Reynolds-averaged Navier-Stokes analysis. A 
D-optimal parametric matrix was selected, and 36 separate CFD 
simulations were included in a full RSM model of RBCC nozzle 
performance, based on 7 sp efficiency. RBCC nozzle performance 
ranged between 77-95% efficient. 

Numerical flow visualization was used extensively to reveal some 
of the underlying physics behind the statistical results. Several inter- 
esting fluid phenomena have been identified for further research. It 
is apparent that the free-expansion process, directly downstream of 
the rocket nozzle, influenced the overall RBCC nozzle performance. 
Maximizing the efficiency of this free expansion will demand some 
creative design solutions. 

A response surface model of RBCC nozzle performance was 
developed, with 13 significant terms affecting the shape and flow 
conditions of a generic RBCC system. Predictions of system effi- 
ciency can now be made throughout our continuous, six-dimensional 


STEFFEN ET AL. 


1039 


design space, without the need for further CFD simulations. Predic- 
tion accuracy to within approximately (±1%) of a similar CFD 
simulation is expected, with 95% confidence. This uncertainty esti- 
mate appears to be valid, in light of the eight additional confirmation 
tests completed. This model can now be used in system design and 
analysis for RBCC propulsion optimization. 

Finally, the cornerstone of the present work was the incorporation 
of CFD into the framework of statistical DOE. The D-optimal ex- 
perimental design provided a parametric matrix of CFD simulations 
that efficiently explored the six-dimensional design space. We assert 
that a comprehensive approach to uncertainty/bias estimation was 
incorporated into the present computational study. We have confined 
our discussion of simulation uncertainty/bias to the issue of verifi- 
cation, or accurate representation of the governing model equations. 
The statistical framework of DOE focused our analysis upon the in- 
terpretation of statistically significant results, using this uncertainty 
estimate as the discriminator. Row parameters integrated from the 
CFD solution (such as mass flow, thrust, and specific impulse) ap- 
pear to be particularly well suited to the regression techniques of 
response surface models. 
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